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Abstract;  Wave  propagation  through  volume  randomness  models  much  ocean  acoustic  behavior 
Random  sound  sp^  introduces  non-linearity  in  the  differential  equations  for  the  ensemble  averages 
This  leads  to  results  which  can  be  significantly  diffeient  than  if  an  averaged  volume  parameters,  such 
as  sound  speed,  are  used  in  a  deterministic  equation.  We  examine  some  theoretical  aspects  of 
randomness  in  the  normal  mode  arJ  parabolic  equatiorts.  As  an  example,  transmission  losses  ir 
several  specific  examples  are  computed  where  differences  result 


INTRODUCTION 


The  study  of  waves  through  random  media  has  been  active  for  about  four  decades  and  has  several 
comprehensive  reviews,  including  Ishimaru  (1978)  and  Ratti  (1983)  (which  emphasizes  ocean 
acoustics).  One  wishes  to  know  about  the  average  solution  to  a  wave  equation  with  urKeruin  sound 
speed.  Frisch  (1968)  points  out  that  this  type  of  problem  is  nonlinear  in  its  dependence  on  the 
coefficients.  As  an  introductory  example,  consider  the  one  dimensional  wave  equation  with  constant 

sound  speed  «„  with  initial  values  «(0.x)  =  /(x).  a,(0.x)*0.  This  problem  has  solutions 

that  can  be  expressed  u(t,x)  =  >^I/(x  -ct)  +  /(x  +  cr)).  A  probability  distribution  of  sound  speeds 

give  averages  of  the  form  (/)  =  |  f  dP{c) »  J  fp(c)dc .  The  average  solution  convolves  initial  values 
with  a  time  dependent  probability  distribution: 


'u)(t.x)=:J/(y)lff,(y)<fy,  Yi(y)-^ 


For  example,  the  normal  distribution  ♦(c)*(crV2>t)"'exp(-(c-Co)V2<T^)  and  the  initial 


(I) 


condition,  /(x)  =  e**,  give  an  average  solution  «(r.x)  =  e‘"  cos4eo/exp<-{Aa/p  /2).  This  shows 
attenuation  of  the  mean  field.  The  average  solution  contains  information  about  the  unceruinty  as 
decay  with  time:  it  behaves  difTerenily  from  the  simple  oscillatory  solution  for  the  averaged  sound 
speed.  One  can  also  show  saturation  phenomenon,  i.e.,  variations  comparable  to  amplitude  as  r 
becomes  large.  In  the  same  case,  but  given  boundary  values  at  x  =0  and  radiation  conditions  for 
X  instead,  one  can  show  phase  speed  of  the  average  solution  is  less  than  the  average  phase 
speed.  Note,  however,  that  tte  solution  is  linear  in  the  initial  conditions,  which  means  the  average  of 
initial  conditions  gives  the  solution  equal  to  the  average  of  solutions  with  those  initial  amdiiions. 


EXACT  SOLUTIONS  TO  RANGE  INDEPENDENT  CASE 


TTw  source  of  the  non-linearity  in  this  case  is  the  inversion  of  the  differential  operator  to  find  the 
solution.  To  see  this,  let  us  consider  the  case  of  a  horizontally  stratified  ocean  model  with  flat 
pressure  release  top  surface  and  bottom  at  -h  which  models  the  trapped  modes  of  the  Pekeris  two  fluid 
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model  (cf  Kinsicr  et  al.  (1982)).  When  this  is  dnvcn  by  a  conunuous  wave,  pouu  source  dus  reduces 
to  a  two-dimensional  Hclmholiz  equation,  which  in  axi-symmeuic,  range  indcpendeni  torm  is 

[r~^d,rd,  +  ))  =  (2?re)"' 5{s  -  >5t/').  (2i 

Boundary  conditions  on  the  top  are  m(i  *  0)*=  0  and  the  bottom  kcosB  uii.-h) -%md  u,{i.~h)  =  0 
This  leads  to  a  completely  discrete  specixum.  Conimuous  specimm  can  be  treated  completely 
analogously,  w  iih  integration  over  the  spcciniro  rather  dun  summauon.  ITiis  assumjHion  precludes 
the  more  complicated  case  where  a  fixed  part  of  the  spectrum  is  discrete  in  some  pans  of  die 
probability  space  and  continuous  in  others.  To  isolate  the  non  linear  step,  use  the  Hankel  iransfunii 
method  in  Ahluwalia  and  Keller  (1977),  Define  the  Hankel  transform  and  its  inverse  by 

H/(j)  =  2»jyo(^i) H’'/(/')  =  -jLjPo(rj)i/(j)dt.  0} 

In  the  Hankel  uansformed  equation  substitute  s  =  ta  to  obtain 

-*•  +a^)j/ =  5(i  -  Zo).  (4) 

The  solutions  may  now  be  expressed  as  inverse  Hankel  transforms  of  the  the  Green  s  furKiion  First 
define  two  fundamental  solutions:  /^(dj  )  which  matches  bottom  boundary  condition  and  /Ja;: ) 

which  matches  top  boundary  condition.  The  Green's  fuiKiion  solution  is 

G(z.Zo)^W-'/,(zJ/^(z^).  (5) 

whercz^  =  min(r,2g).  z,  =  max(z,rg)  and  the  Wronskian  is  W  =  W(tn)  = /, , /j -/, /j , .  Thus  the 
solution  to  wave  equation  is  the  inverse  Hankel  transform  of  Green’s  funaion 

u(r,z)  =  H"'G.  (6) 

(6)  uses  a  single  z-dependem  sound  speed  profile.  The  average  of  solutions  for  many  ptofilcs  is 

(u)  =  £(H-'C)  =  H-'£(C)  =  ir'{w-'/,  /j).  (7) 

£  (C )  is  an  averaged  Green’s  function,  and,  obviously,  one  cannot  simply  pull  the  averages  into  its 
numerator  and  denominator  separately.  One  could  deform  contour  of  inverse  Hankel  transform  to 
obtain  ’branch"  cuts  instead  of  poles  for  the  normal  mode  contributions.  However,  the  same  result 
follows  from  normal  mode  theory  with  averaging. 

NORMAL  MODE  REPRESENTATION 

Ahluwalia  and  Keller  (1977)  obtain  the  range  independent  rxirmal  mode  representation  of  solution 

=  (8) 

It  has  terms  satisfying  the  right  radiation  condition  if  the  radial  function  is  ^*(r)=  The 

coefficients  depend  on  source  depth  actxrrding  to  =  ^mCzo)  where  vertical  equation  is 

l^f  +  k^(n^+aj)}^^(z)  =  0.  (9) 

The  asymptotic  behavior  of  vertical  eigenvalues  is 

--f|■^(2xrl^  +  0)^  m-4-.  (10) 

Set  q-k^n^.  Qearly  the  eigenvalues  a*  =  a*  will  change  as  q  changes.  One  can  quantify  this  as  a 
functional  derivative.  Fust,  rewrite  the  z  equation  (9)  ( '  is  r  difTerentiation) 

0;+(?-a,)^,=O  (11) 

Let  q  vary  according  \oq=9q  +  ep  ,  then  one  can  compute  the  Frtehet  derivative  (see  Ito  et  al. 
(1987)) 


Kr  =  P)  =  <^^(<7)1  P 1  = 


«?£ 


(12) 


aoo 
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DifTcreniiatc  entire  equation  (11)  with  respect  to  e  to  obtain 

^  +  (P-  =t>  ‘  i 

.0 

Tfic  inner  product  is  (  f,g)  =  j  ^f(z)g(z)di  .  For  iwmialized  eigenvalues,  one  dillerenuaics  il*r 
normality  condition  to  find  that 

0  =  .  (14, 

Take  inner  pnoduct  oi  with  diffcrendaicd  equation  { 1 3).  and  intcgralc  by  pans  with  respect  to  ; 

Since  the  boundary  conditions  for  are  constant,  they  carry  over  to  This  gives  Fttcltci 
derivative  of  eigenvalues 

da(q){p\-a^^=(p^^,p^).  (15) 

As  a  simple  demonstration,  take  sound  speed  constant  <7  =  0,  6{0)  =  6{k)  =  0 ,  k  =  1 .  and  h~x  , 

then  a  =  =  -1  and  6 -mi.  (15)  becomes  da^(qip\  =  2  /  )r|*p(r)sin^m2di. 

Also,  one  may  apply  (15)  to  compute  average  contribution  from  each  mode,  if  the  probability 
distribution  of  the  speeds  can  be  reparameirized  by  a  =  a„ .  U  q-q(<o).  p  =  3q  /  Sw  .  then  one 
may  compute  the  a  dependence  by  dct*(flXpIdlP(a,)  =  dP(o)). 

Turn  now  to  the  calculation  of  the  average  mode's  contribution  when  the  sound  speed  is  constant,  but 
statistically  uncertain.  This  contribution  will  now  depend  on  many  horizontal  wave  numbers  rather 

than  just  one.  If  6,  =(/j«  +  e)//i  and  a,  =  •yk^-(fljr  + 0)*/i'^  ,  the  series  (8)  becomes 

sinb,i  W^'^o,r),  (16) 

For  a  rigid  boaom  condition  pick  0s=l/2.  Take  a  uniform  probability  distribution  over  a  smaU 
interval:  dP  (k  )=l/(2c).  ko-£<k<kQ  +  c.  Three  cases  arise  according  10  vvhether  the  eigenvalue 
for  the  probabilistically  distributed  mode  is  always  positive,  always  negative,  or  spread  around  zero. 
The  last  case  includes  the  previous  two,  so  consider  it  Suppose  that  the  sound  spetd  distribution  is 

such  that  the  eigenvalues  are  in  a  small  band  centered  about  zero:  k^nx*0  )lh  ,  k  =  k  -  k^.  One 
finds  o,  -  [2k<,»:]‘'^^.  Suppose  r  is  large  enough  that  one  may  use  the  asymptotic  form 

For  a  single  term  of  (16),  only  the  Hankei  function  is  affected  by  variation  in  eigenvalues.  Thus  the 
average  of  a  term  of  ( 16)  can  be  written  as  a  product  of  the  sines  times  an  integral  involving  the 

Hankei  function.  Let  s^r^lk^K,  that  integral  becomes 

I  ,1,}  = 

These  incomplete  gamma  functions  appear  instead  of  Hankei  functions  for  yr,  in  (16). 

Uncertainty  in  the  bottom  condition. in  the  equations  (16)  give  another  example.  One  may  model  the 
uncertainty  in  the  difference  between  the  velocity  of  two  fluid  layers  by  random  boaom  boundary 
condition.  S&.k-OtK+S^)/h,dP($)=ddfZe,  6Q-£<e<6o+£,  ht=:6-6o.  Thus*,  =/t  +  /. 


